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摘要 :以 山西 太 岳 山 华 北 落叶 松林 地 为 主 的 小 流域 作为 研究 对 象 ,采用 地 统计 学 方法 结合 地 理 信息 系 统 (GIS) 技 术 手 段 ,研究 
9 的 空间 变异 特征 ,以 及 植被 分 布 和 地 形 因子 对 其 影响 规律 @ 结果 表明 ;在 时 间 稳 定性 的 前 
提 下 ,土壤 水 分 含量 和 变异 系数 随 土屋 加 深 逐 渐 降低 。 三 层 土 壤 水 分 半 方 差 函数 的 最 优 拟 合 模 型 为 球状 模型 , 变 程 范围 在 
1.1 一 1.4 km 和 的 空间 自 相关 性 ,其 中 0 一 20 cm 和 20 一 40 cm 层 土壤 水 分 的 室 幅 虽 质 性 程度 高 于 40 一 60 cm 土 层 ,以 
中 间 层 的 结构 因素 占 总 变异 比例 最 大 。 自 然 结 构 因 素 ( 地 形 . 母 质 .植被 和 土壤 等 ) 对 不 同 土 层 土壤 水 分 的 总 空间 变异 性 起 主 
导 作 用 (81.4% 一 91.3% ) ,而 随机 因素 (取样 误差 、 人 为 干扰 等 ) 的 影响 相对 较 小 (8.7% 一 18.6%)。 沿 着 集 水 线 由 西 -东方 向 ,从 
边缘 的 土壤 水 分 高 值 斑 块 区 逐渐 过 渡 到 明显 的 低 值 斑 块 区 ,梯度 变化 明显 J 人 研究 发 现 ,在 植被 覆盖 异 质 性 小 的 山地 ,土壤 水 分 
的 空间 蜡 质 性 主要 由 地 形 因 素 引 起 ,具体 表现 为 其 与 坡 向 指数 (TRASPD》 、 城 度 海拔 和 土壤 有 机 磋 、 全 氮 呈 极 显 著 相关 关系 
(P < 0.01) ,而 与 植被 指数 (NDVI) 呈 弱 的 负 相 关 关 系 。 闪 加 分 析 显 示 ) 在 阴 坡 、 坡 度 较 缓 ( < 15°) 及 高 海拔 车 合 的 区 域 土壤 水 
分 含量 较 高 。 研 究 结 果 可 为 山地 人 工 林 构 建 和 植被 恢复 中 土壤 水 资源 的 利用 以 及 水 分 管理 策略 的 制定 提供 理论 依据 。 
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Abstract: Soil moisture is a hydrological characteristic that plays a crucial role in the soil-plant-atmosphere continuum. The 
spatial distribution of soil moisture and factors that influence this distribution are extremely important to many ecological 
processes, especially in temperate mountain forest regions such as those in northern China. Nonetheless, little is known 
about soil moisture heterogeneity in typical subalpine forest lands or the controlling factors of soil moisture change. In this 
study, sOil moisturescontent (SMC) was measured for depths of 0 to 60 cm in a typical watershed forest, dominated by 
Larix principis-ruppreehtii, in Taiyue Mountain, Shanxi Province. We used Geographic Information System (GIS) data and 
géostatistiéal methods to examine the spatial heterogeneity of SMC and the factors (i.e., vegetation and topography ) that 
inifluence it. SMC and its coefficient of variation ( CV) decreased with increasing soil layer depth, exhibiting high to 
moderate variation under an assumption of time stability. At this small watershed scale, the best-fit model for SMC of the 
three layers (0 一 20 cm, 20—40 cm, and 40 一 60 cm) is the spherical model; the spatial autocorrelation distance ranged 


from 1.1 to 1.4 km. The distribution of SMC has strong spatial autocorrelation in all three layers. The degree of spatial 
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heterogeneity of soil moisture in the 0 一 20 cm and 20 一 40 cm soil layers was higher than that in 40—60 cm layer. The 
structural characteristics of the middle layer accounted for the largest proportion of the total variation. Spatial variability of 
soil moisture was mainly (81.4% to 91.3% ) controlled by natural structure factors such as topography, parent material, 
vegetation properties, and soil properties, while the influences of random factors such as sampling error, human 
disturbance ，plant root distribution, and absorption were relatively small (8.7% to 18.6% ). Moran's J values indicated that 
soil moisture had moderately aggregated distributions ( 0.44 一 0.51 ) across soil depths. The overall trend of SMC spatial 
distribution gradually transitioned from high to low values from west to east along with the waterline, an obvious gradient 
change. We also found that the spatial heterogeneity of soil moisture was mainly controlled by terrain factors in this 
mountainous catchment with homogeneous vegetation cover. Spatial variability of soil moisture was significantly correlated 
with topographic factors (such as aspect, slope, and elevation ) and soil chemical properties (such as soil/organic carbon 
and total nitrogen) (P<0.01). However, it was weakly and negatively correlated with the normalizéd differéenée vegetation 
index (NDVI). Overlay analysis suggested that SMC was relatively high in regions where the thrée attiibutés of shady slope， 
smooth slope (<15 degrees ) ，and high altitude overlapped. Moreover, soil moisture variation néarly matched the color 
depth change in a superimposed terrain map. Taken together, the results of this study provide data,for planning the use of 
soil water resources in artificial vegetation restoration and construction of temperateAubalpine forest and a theoretical basis 


for devising water management strategies. 


Key Words: small watershed ; soil moisture content; spatial heterogeneity; Larix principis-rupprechtii; influencing factors 


水 是 维持 植物 生长 的 主导 因子 ,是 维持 森林 生态 系统 稳定 性 的 重要 因素 之 一 '" 。 在 流域 尺度 上 ,水流 在 
重力 作用 下 汇集 ,使 得 不 同 景观 要 素 之 间 的 物质 流 和 能 量 流 的 联系 更 为 紧密 ,而 流域 的 下 热 面 因素 (如 地 形 、 
植被 和 土壤 特征 等 ) 均 会 对 流域 的 水 文 特征 产生 最 著 的 影响 。 在 暖 温带 山地 针 叶 林 区 ,土壤 水 分 作为 流域 
水 文 特征 的 关键 表征 ,是 森林 生态 系统 物质 循环 的 载体 ,在 土壤 -植被 -大 气 连续 体 间 水 文 过 程 及 能 量 交换 中 
起 着 重要 作用 '。 由 于 土壤 水 分 与 地 形 <《 土 壤 和 和 植被 间 的 复杂 作用 关系 ,其 在 不 同 空 间 尺 度 上 存在 较 大 变 
异性 。 

近年 来 ,国内 外 学 者 对 于 土壤 水 分 不 同 尺 度 的 空间 异 质 性 进行 了 大 量 研究 ,地 理 信息 系统 和 地 统计 学 方 
法 相 结合 已 成 为 此 类 研究 的 有 效 工 具 , 但 主要 集中 在 农田 和 草地 ,有 关于 林地 土壤 水 分 空间 变异 的 研究 较 
少 ! 。Savva 等 站 研究 \ 了 美国 马里 兰州 郊区 不 同 森林 和 草地 覆盖 下 浅 层 土壤 水 分 的 空间 分 布 格局 。 
Rosenbaum 等 5 在 德国 奴 杉 覆盖 为 主 的 集 水 区 展开 研究 ,揭示 了 不 同 深度 (5.20 .50 cm) 土壤 水 分 的 时 空 变异 
规律 。 王 存 国 等 对 长 白山 阔 叶 红 松林 3 块 50 mx50 m 样 地 表层 土壤 水 分 空间 异 质 性 进行 了 研究 。 上 述 研 
究 的 范围 从 样 地 到 集 水 区 尺度 不 等 ,以 低 海拔 (<800 m) 的 土壤 表层 水 分 研究 为 主 。 对 于 亚 高 山 山地 森林 覆 
盖 为 主 的 / 屏 流 域 人 研究 甚 少 ,地 形 的 起 伏 使 得 土壤 水 分 的 再 分 配 过 程 更 加 复杂 ,而 土壤 更 深层 次 水 分 的 空间 变 
异 对 于 植被 的 生长 和 生态 系统 水 分 平衡 同样 具有 重要 的 作用 。 

国内 在 平 早 半 干 旱 黄土 丘陵 区 的 小 流域 研究 发 现 , 土 壤 水 分 空间 异 质 性 是 多 重 太 度 上 环境 因子 共同 作用 
的 结果 ,降水 强度 地质 地貌 等 主要 影响 其 空间 分 布 “ 中 。 国 外 研究 表明 ,小 流域 土壤 水 分 受到 土壤 性 质地 
形 特征 .气象 条 件 .植被 覆盖 和 地 下 水 状况 等 因素 的 影响 '" ,但 对 于 不 同 的 空间 尺度 ,各 环境 因子 对 于 其 影响 
程度 不 尽 相同 ,景观 尺度 上 土壤 水 分 的 时 空 异 质 性 与 环境 因子 的 关系 已 成 为 国内 外 的 研究 热点 号 。 华 北 晋 
东南 的 亚 高 山 山 地 主要 以 针 叶 人 工 林 为 主 , 有 学 者 提出 大 规模 的 造林 虽 增 加 了 碳 储量 ,但 可 能 会 引起 土壤 干 
燥 化 效应 2。 因此 ,掌握 土壤 水 分 的 空间 尺度 分 异 特征 及 其 影响 因素 对 于 华北 森林 景观 的 有 效 经 营 与 管理 ， 
防止 土壤 干 化 与 侵蚀 ,实现 森林 生态 服务 功能 的 最 大 化 发 挥 都 显得 尤为 重要 。 本 研究 从 小 流域 扩 度 上 对 暖 温 
带 森 林 土 壤 水 分 的 空间 分 布 模式 、 尺 度 依赖 性 及 其 控制 因子 进行 研究 ,为 理解 亚 高 山 山地 森林 不 同 尺 度 上 土 
壤 水 分 的 空间 变异 及 水 文 过 程 提供 依据 。 
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1 研究 区 与 研究 方法 


1.1 研究 区 概况 

研究 地 处 于 山西 省 太 岳 山 好 地 方 林场 的 接 石 沟 集 水 区 (36°47'24 一 43”"N ,111°59'43” 一 112°03'34"E) , 面 
积 为 5.65 km ,隶属 于 太行 山 支脉 太 岳山 中 部 ,黄土 高 原 的 东 缘 。 该 地 区 气候 属 暧 温带 半 湿 润 大 陆 性 季风 气 
候 , 年 平均 气温 8.7 %C ,年 日 照 2600 bh, 年 平均 降水 量 为 660 mm ,雨量 主要 集中 在 7、8.9 月 ,无霜期 110 一 480 
d。 集 水 区 海拔 2037 一 2359 m ,岩石 主要 为 石灰 岩 和 花 岗 片 麻 岩 , 土 壤 亚 类 有 褐 土 、 标 壤 和 山地 草 名 区 随 海 
拔 呈 地 带 性 分 布 。 

研究 区 内 ,植被 类 型 以 华北 落叶 松 ( Zarix principis-rupprechtii) 人 工 林 为 主 , 约 占 流域 面积 的 80.5% ,其 中 以 
近 熟 林 (31 一 40a) 为 主 , 广 泛 分 布 于 集 水 线 两 侧 。 次 生 林 如 白桦 ( Berula plaiyphylla) 和 黄 华 柳 @S4lixheuprea) 等 
零星 分 布 于 小 流域 海拔 较 低 的 区 域 , 常 与 年 龄 较 小 的 落叶 松 (<30a ) 形 成 针 阔 混 驳 林 。 各 植被 尖 型 的 林 分 结 
构 特 征 见 表 1。 灌 木 主 要 有 土 庄 绣 线 菊 (Spiraea pubescens ) 、 毛 栋 ( Corylus mandshiriea ) 骨 枝子 (Lespedeza 
jicolor) 、 黄 刺 玫 ( Rosa xanipina) 等 ;草本 层 的 代表 植物 有 细 叶 蔓草 ( Carex Ilancedlat@) 、 小花 风 毛 菊 ( Saussurea 
parviflora) ,小 红 菊 ( Dendranthema chanetii) 等 。 


表 1 研究 区 主要 植被 类 型 的 林 分 结构 特征 


Table 1 Stand structure characteristics of the main vegetation types in the study area 


密度 平均 树 高 平均 胸径 平均 宁 旺 平均 冠 高 闭 度 
林 分 类 型 林 龄 设 均 树 高 均 胸径 均 胶 幅 均 冠 高 郁 闭 度 土 层 深度 
Qaridl AZ Density / Mean Mean Mean crown Mean crown Canopy Soil depthy 
ee Se ( 株 /hm?) height/m DBH/Acm width /Am height /m density/ % eh 
华北 落叶 松 人 工 林 ” 幼 龄 林 ( <20) 1566.67 8.19 58 3.29 7.21 80.00 0.50 
中 龄 林 (21 一 30) 1027.78 13.09 13.40 4.34 10.10 75.00 0.55 
近 熟 林 (31 一 40) 1576.73 16.47 20.67 4.95 13.60 80.00 0.60 
成 熟 林 ( >40) 1238.89 20.25 25.87 5.95 16.10 85.00 0.75 
阔 叶 林 一 1677.83 7.15 7.38 3.14 5.63 80.00 0.50 
针 阔 混交 林 = 2016.67 6.38 6.98 2.73 4.15 85.00 0.55 


1.2 样品 采集 与 测定 

结合 研究 区 1:5 万 地 形 图 和 和 植被 分 类 图 ,于 2015 年 10 月 20 日 采集 研究 区 土壤 样品 。 因 高 海拔 山区 地 形 
复杂 .采样 难度 较 大 ,常用 的 网 格 采样 法 在 本 研究 区 实施 困难 ,因此 采用 代表 不 同 植被 类 型 .地形 条 件 ,兼顾 均 
勺 分布 的 采样 方法 选择 样 地 中 ( 表 2)。 样 点 数量 为 81 个 ,采样 间隔 约 为 260 m, 确 保 均 匀 取 样 ,充分 涵盖 小 
流域 的 自然 特征 。 海 拔 . 堪 度 . 坡 向 等 地 形 因子 从 DEM (Digital elevation model) 中 提取 ,DEM 由 10m 等 高 距 的 
地 形 图 数字 化 产生 。 植 被 类 型 分 布 从 2014 年 9 月 30 日 高 分 辨 率 的 资源 三 号 遥感 卫星 影像 解 译 获 取 。 另 外 ， 
植被 生长 状态 和 楠 被 覆盖 度 用 归 一 化 植被 指数 NDVI 表示 ,该 指数 从 影像 多 光谱 数据 的 近 红 外 波段 和 红外 波 
段 计算 而 得 弟 

为 了 便于 分 析 比 较 , 在 取样 前 一 周 内 天 气 稳定 ,基本 无 明显 降雨 的 晴天 集中 采样 ,并 于 1 天 内 尽量 将 土 样 
采集 完毕 。 用 手持 GPS 记录 采样 中 心 点 坐标 ,在 以 该 点 为 中 心 的 Im 正方 形 小 样 方 内 均匀 确定 3 个 采样 点 ， 
每 淮 采 样 点 用 直径 为 3 cm 的 土 钻 取 表层 (0 一 20 cm) .中层 (20 一 40 em) 底层 (40 一 60 em) 共 3 层 土壤 ,每 层 
分 别 取 20 g 土 样 , 噜 除 明 显 的 植物 断根 和 枯 落 物 等 人 杂质 ,等 量 充分 混合 后 取 1/3 装 入 编号 的 馈 盒 中 , 带 回 实 
验 室 。 采 用 同样 的 方法 另 取 100 g 土 样 装 入 自封 袋 测定 土壤 的 化 学 性 质 。 表 层 土 共 采集 162 份 样品 ,但 由 于 
部 分 采样 点 的 基 岩 深度 较 浅 ,中 层 和 底层 的 土 样 数量 分 别 为 144 份 和 120 份 。 采 用 烘 干 法 对 土 样 的 土壤 含水 
量 进行 测定 ,将 称 量 后 的 土壤 样品 (精确 到 0.01 g) ,在 (105+2)% 的 烘箱 内 烘 至 12 h, 取 出 后 称 量 ,进而 计算 
土壤 含水 量 , 共 测定 81 个 样 点 的 213 份 样品 。 另 外 同等 数量 的 样品 采用 Flash 2000 有 机 元 素 分 析 仪 (Thermo 
Scientific ，Brookfield ，Wisconsin，USA ) 测定 土壤 有 机 碳 和 全 氮 。 
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表 2 研究 区 不 同 地 形 与 植被 类 型 取样 点 数量 分 布 


Table 2 ” Distribution of sampling points in different terrain and vegetation types 


分 类 植被 类 型 样 点 数 海拔 样 点 数 坡 向 样 点 数 坡度 样 点 数 
Class Vegetation types Number Elevation A/m Number Aspect Number Slope/(°) Number 
1 华北 落叶 松 幼 龄 林 13 <2200 39 阴 坡 19 <10 4 
2 华北 落叶 松 中 龄 林 20 2200 一 2300 14 半 阴 坡 26 10 一 15 11 
3 华北 落叶 松 近 熟 林 27 >2300 28 阳 坡 26 15 一 20 24 
4 华北 落叶 松 成 熟 林 8 半 阳 坡 10 20 一 25 20 
3 阔 叶 林 3 >25 22 
6 针 冰 混交 林 6 
7 灌木 林 4 


”华北 落叶 松 幼 龄 .中 龄 . 近 熟 和 成 熟 林 分 别 是 <20a, 21 一 30a，31 一 40a 和 >41a 生 的 华北 落叶 松 人 工 林 , 阔 叶 林 主要 为 让 梯次 歼 林 , 针 间 混交 
林 为 华北 落叶 松 和 白桦 混交 林 
在 数据 处 理 时 , 坡 向 方位 角 按照 Roberts & Cooper :的 方法 进行 了 转换 , 坡 向 由 0 会 360* 的 罗盘 值 转换 成 
一。 0 一 1 之 间 的 值 。 转 换 公 式 为 : 
TRASP={1-cos[(m/180) (aspect-30)1}/2 

式 中 ,TRASP (transformation of aspect) 为 坡 癌 指数 ;aspect 为 坡 向 方位 角度 。 通 过 转换 ,TRASP 介 于 0 一 1 之 
间 , 其 值 越 大 ,表示 生境 越 干 热 ,0 代表 北 北 东方 向 ,1 代表 南 南西 方向 。 
己 。 13 空间 格局 分 析 
( ”1.3.1 半 方 差分 析 
OO 林地 土壤 水 分 为 区 域 化 变量 ,同时 具有 随机 性 和 结构 性 《空间 自 相 关 性 ) 特征 , 当 土 壤 水 分 变量 满足 二 阶 
0 平稳 假设 和 本 征 假设 , 且 样 本 数量 足够 大 时 ,可 应 用 半 方 差分 析 5 。 根 据 Matheron™"] , 半 方 差 函 数 y(h) 计 
中 。 算 公式 为 : 


NOQA) 

PT [z(x,) — z(x, +h)] 
a 式 中 ,hh 为 两 样本 点 空间 分 隔 距 离 ; N(h) 为 在 分 隔 距 离 为 h 时 的 样本 对 总 数 ; z(x;) 和 z(x; + 及) 分 别 为 区 域 
9 化 变量 z(x) 在 空间 位 置 x 和 x + 全 处 的 实测 值 。 
一 半 方 差分 析 结 果 反 映 了 数据 的 空间 关联 程度 ,只 有 数据 具有 空间 自 相 关 性 ,空间 插值 才 有 意义 。y(h) 中 
) ”的 块 金 值 (C,) 表 示 随 机 部 分 的 案 间 异 质 性 ,而 (C) 表 示 空 间 自 相关 部 分 的 空间 异 质 性 , 基 台 值 (C6+C) 表示 区 
域 化 变量 的 最 大 变 虹 。 屎 CwGGN+C) ] 为 基底 效应 ,表示 样本 间 随 机 因素 引起 的 变异 程度 ,该 值 越 小 ,空间 相 
关 性 越 强 ,按照 区 域 化 变量 空间 相关 性 程度 的 分 级 标准 ,该 值 为 <25% .25% 一 75% 、>75% 时 ,分 别 表明 变 
量 的 空间 相关 性 为 较 强 \ 中 等 . 较 弱 。 
1.3.2 空间 自 相关 分 析 

空间 格局 的 关联 程度 采用 了 Moran 的 了 指数 法 ,采用 该 指数 以 探 明 土壤 水 分 变量 的 分 布 是 否 聚 集 .离散 
或 者 随机 分 而 。 公 式 为 [5 : 


y(h) = 


Moran’s 1 = 一 一 
人 Er 3 > (x — #)? 
式 中 ,x; 和 x 分 别 为 景观 要 素 x 在 空间 单元 i 和 j 中 的 观测 值 ， kg 的 平均 值 , WW; 为 相应 权重 ,N 为 空间 单元 
总 数 。Moran's 7 的 值 域 为 [ -1,1] , 取 值 为 -1 表示 完全 负 相 关 , 取 值 为 1 表明 完全 正 相关 ,而 取 值 为 0 表示 随 
机 分 布 ,无 空间 自 相 关 性 
1.4 数据 处 理 与 分 析 

运用 SPSS18.0 对 数据 进行 Kolmogorov-Smironov (K-S ) 正 态 分 布 检验 和 基本 统计 分 析 ,运用 地 统计 学 软件 
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包 GS+ for windows 9.0 和 Arcgis10.0 进行 半 方 差分 析 和 Kringing 内 插 局 部 估计 。 


2 结果 与 分 析 有 mm 降雨 量 一 气温 一 一 取样 时 间 点 


[ea 
un 


2.1 土壤 水 分 的 取样 时 间 与 统计 特征 


研究 区 自动 气象 站 2015 年 的 降雨 监测 数据 显示 | 3 
(图 1),10 月 份 降雨 总 量 为 58.67 mm, 月 均 温 为 “号 1 
5.08%C 。 降 雨 主要 集中 在 月 底 ,本 次 采样 的 时 间 点 为 10 时， 
月 20 日 ,在 18 日 只 有 一 次 0.25 mm 的 小 雨 ,加 之 气温 站 
较 高 , 蒜 发 量 大 ,土壤 较 干 ,为 分 析 土 壤 水 分 的 空间 异 质 “| 
性 提供 了 时 间 上 的 相对 稳定 性 。 表 3 为 小 流域 各 植被 " 2015-10.05 2015-10-12 2015-10.17 2015-10-26 


类 型 下 不 同 土 层 中 土壤 水 分 的 描述 性 统计 结果 。 从 表 日 期 Date 
中 可 知 , 随 着 土 层 深度 的 逐渐 加 深 , 土壤 含水 量 的 平均 1 研究 区 10 月 信 m 攻 这 特 贡 到 栏 时 间 上 
值 依 次 降低 。 不 同 土 层 的 土壤 水 分 的 水 平 变异 均 表现 Fig.1 Characteristics of rain and heat of October in the study 
为 自 表 层 向 下 显著 降低 趋势 ,土屋 40 cm 以 上 表现 出 强 。 area and the sampling Eapointg 
变异 性 (>35%),40 cm 以 下 为 中 等 变异 性 (15% < 
CV<35%)'™ 

半 方 差 廿 效 通 过 寸 区域 化 变量 分 割 等 距离 样 点 间 的 差异 来 研究 变量 的 空间 相关 性 空间 相关 性 分 析 的 
变量 须 满足 正 态 分 布 !'5 。 表 3 中 显示 K-S 检验 的 he te de 
析 , 并 且 20 em 以 下 土 层 土壤 水 分 含量 的 已. 值 要 明显 优生 0 全 20 cm 土 层 。 


表 3 不 同 土 层 深度 土壤 水 分 的 描述 性 统计 结果 


Table 3 Classical statistics analyses of the soil moisture in different soil layers 


举 本 量 本 示 准 差 恋 异 系 间 同 
十 层 “不 屋 最 小 信 最 大 信 AAA 人 SEE 关 。 3 夺 系 ， 太朗 峰 度 KES 值 
Depth/ cm 一 Min./% Max./% Mean XY % ey es Skewness Kurtosis K-S value 
size deviation coefficient /% 
0 一 20 81 5.46 28.90 15.75 6.14 38.98 0.34 -0.73 0.29 
20 一 40 72 4.89 25.76 13.02 4.81 36.94 0.35 0.47 0.85 
40 一 60 60 5.66 21.67 12.54 4.09 32.62 0.45 0.64 0.60 


不 同 土 层 样本 量 不 一 致 是 由 各 梯 点 忆 层 这 度 不 同 所 致 


2.2 不 同 深度 土壤 水 分 的 空间 异 质 性 

图 2 和 表 4 为 不 同 士 层 深度 土壤 水 分 在 各 向 同性 下 的 变异 函数 理论 模型 及 相关 参数 。 从 表 中 可 知 ,不 同 
ee 天国 奉公 放 册 上书 用 的 3 个 层次 均 具 有 明显 的 空间 自 相关 结构 ， 
雪 P0 一 20cm 和 20 一 40cm 层 土壤 水 分 的 基 台 值 高 于 40 一 60cm 土 层 , 说 明 表 层 和 中 层 的 空间 蜡 质 性 程度 


表 4 各 向 同性 下 变异 函数 理论 模型 及 参数 


Table 4 Parameters of the best-fitted semi-variogram model for isotropic variogram 


决定 系数 块 金 / 基 台 交叉 验证 
论 模 型 全 基 台 恋 程 > 关 a 
竹 层 理论 模 天 金 什 值 变 程 Determination 残 差 Nugget/ Sill Cross-validation 葛 兰 指数 
Theoretical Nugget Sl Range 和。 Residual 7 
Depth/ em coefficient Co/ (Cot Moran J 
model Co CotC A/m 人 RSS RMSE AME 
R’ C) /% 
0 一 20” Spherical 0.031 0.190 1140 0.90 4.40x10™3 16.5 0.83 0.023 0.492 
20 一 40” Spherical 0.016 0.183 1393 0.83 8.32x10-3 8.7 1.11 -0.011 0.436 
40 一 60” Spherical 0.020 0.110 1115 0.68 5.35x10™3 18.6 1.03 0.027 0.511 


* 表示 数据 经 log 转换 ,因数 据 经 log 转换 后 ,空间 结构 的 规律 性 更 明显 , 量 残 差 较 小 ,符合 正 态 分 布 ;RMSE 为 标准 均 方 根 预测 误差 ( Root-Mean-Square 
Standardized Error) , 值 越 接近 于 1 ,误差 越 小 ; AME 为 平均 标准 误差 ( Average Mean Error) , 值 越 接近 于 0, 误差 越 小 
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Ea @D 

术 .S 

+ 

0 531.39 1062.79 1594.18 
口 口 
玫 吕 = 
图 呈 口 
糙 日 
入 己 口 口 
站 电 
20 一 40cm 40 一 60cm 
531.39 1062.79 1594.18 0 531.39 1062.79 1594.18 


步 长 (h) 


Separation distance/m 


图 2 森林 土壤 水 分 的 各 向 同性 半 方 差 函数 


Fig.2 Isotropic semi-variogram of forest soil moisture content 


明显 高 于 底层 。 三 层 土壤 水 分 的 半 方 差 图 数 与 步 长 尹 的 关系 均 符合 球状 模型 ,具体 特征 表现 为 空间 相关 性 随 
距离 的 增加 而 衰减 ,其 空间 结构 是 当 样 点 间隔 距离 达到 变 程 (如 4。yo .6140 m) 之 前 , 样 点 的 空间 相关 性 随 
样 点 距离 的 增加 逐渐 降低 直至 消失 (图 2)。 

从 决定 系数 RR 和 残 差 RSS 来 看 ( 表 4) ,0 一 20 cm 和 20 二 40 cm 和 层 土壤 水 分 半 方 差 函 数 与 步 长 hh 关系 
的 理论 模型 解释 效率 相对 较 高 (90% 和 83% ) ,而 40 一 60sem 生 层 的 半 方 差 函 数 模型 解释 效率 较 低 (68%),3 
个 土 层 水 分 模拟 的 残 差 RSS 值 均 较 低 (4.40x10” 一 8.32X10 了 7) 。 

由 表 4 可 知 ,不 同 土 层 土壤 水 分 的 空间 异 质 性 中 自 相 关 和 随机 部 分 所 占 比 例 存 在 一 些 差异 。0 一 20 cm 
和 40 一 60 cm 土 层 土壤 水 分 的 Cu,A( Cu+C) 值 相对 于 20 一 40 cm 层 较 高 ,说 明 表 层 和 底层 由 随机 部 分 引起 的 空 
间 异 质 性 高 于 中 层 , 即 中 层 由 空间 自 相关 部 分 引起 的 空间 异 质 性 较 高 。 结 构 方 差 比 C/A( C6+C) 表 明 结 构 因 素 
对 于 表 中、 底层 的 空间 变异 性 的 贡献 率 分 别 达 到 83.5% ,91.3% 和 81.4% , 上 且 残 差 RSS 较 小 。 根 据 相 关 性 程度 
的 分 级 标准 ,三 层 土壤 的 CuX(Co+GE9i 值 均 明 显 低 于 25% ,说 明 具 有 强烈 的 空间 相关 性 。 结 合 当 地 华北 落叶 松 
及 生境 的 实际 情况 , 自然 结构 因素 (如 气候 植被 . 母 质 .地 形 和 土壤 类 型 等 ) 对 不 同 土 层 土壤 水 分 的 总 空间 变 

异性 起 主导 作用 (81.4% 一 -94.3%) ,而 随机 因素 (如 取样 误差 ,人 为 干扰 ,植物 根系 分 布 与 吸收 等 ) 的 影响 相对 
较 小 (8.7% 一 18.69% 。 

变 程 4 可 反映 区 域 化 变量 的 空间 作用 范围 ,其 大 小 受 观 测 尺度 限定 ,在 变 程 范围 内 , 样 点 间 的 距离 越 小 ， 
其 空间 相关 性 越 大 。 由 图 2 和 表 4 可 知 ,0 一 20 cm 和 40 一 60 cm 土 层 土壤 水 分 的 变 程 比较 接近 (1140 m 和 
1115 m) ,而 20 一 40 cm 中 层 土壤 水 分 的 变 程 较 大 (1393 m) ,说 明 其 在 较 长 距离 内 存在 空间 结构 异 质 性 ,而 当 
超过 此 变 程 时 , 注 近 水 分 的 空间 相关 性 消失 。Moran' 1 指数 表示 全 局 的 空间 自 相关 程度 ,3 个 层次 土壤 水 分 的 
Moran 对 (0.436 一 0.511) 均 远 远 大 于 0, 说 明 全 局 土壤 水 分 的 空间 相关 性 较 明 显 , 呈 聚集 分 布 。 其 中 ,底层 土 土 
壤 水 汾 的 Moran’ 7 指数 为 0.511 ,空间 相关 性 又 较 表 层 和 中 层 略 强 。 

2.3 土壤 水 分 的 空间 分 布 格局 及 驱动 因子 分 析 

图 3 为 不 同 土 层 土壤 水 分 的 空间 分 布 特征 。 据 图 可 知 , 沿 着 集 水 线 由 西 -东方 向 ,0 一 20 cm 和 40 一 60 cm 
层 从 边缘 的 土壤 水 分 高 值 斑 块 区 ,逐渐 形成 明显 的 低 值 斑 块 区 ,梯度 变化 明显 。 而 20 一 40 em 的 水 分 高 值 区 
向 下 往 南 部 延伸 ,但 总 体 趋势 与 上 下 两 层 相似 ,呈现 出 西 侧 大 于 东 侧 的 现象 。 结 合 地 形 因 子 ,在 阴 坡 、 坡 度 较 
缓 (<15" ) 及 高 海拔 受 合 的 区 域 土壤 水 分 含量 较 高 ,地 形 因子 到 加 的 高 低 区 域 与 表层 土壤 水 分 的 变化 规律 大 
致 吻合 。 其 中 ,表层 十 最 高 值 区 的 土壤 含水 量 为 24.66% 一 28.90% ,最 低 值 区 为 5.47% 一 7.80% ; 中 层 十 最 高 值 
区 含水 量 为 20.93% 一 25.77% ,最 低 值 区 为 4.89% 一 6.86% ;底层 土 高 值 区 的 含水 量 为 19.07% 一 21.68% , 低 值 
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区 为 5.66% 一 7.27% ,表现 出 全 局 土壤 水 分 的 高 值 区 随 士 层 加 深 逐 渐 降 低 ,而 低 值 区 却 呈 现 出 先 降低 后 升 高 的 


0 一 20cm 20 一 40cm 

土壤 水 分 含量 /% 土壤 水 分 含量 /% 
Ml 5.466—7.795 4.893—6.859 
7.795—9.702 EE 0.859 一 8.316 
9.702—11.26 ul 8.316—9.396 
El 11.26 一 12.54 Wm 9.396 一 10.20 
MM 12.54 一 14.10 本 一 10.20—11.27 
14.10—16.01 Ey 11.27—12.73 
mn 16.01 一 18.34 Wm 12.73—14.70 
El 18.34—21.18 my 14.70—17.35 
Bl 21.18—24.66 Ey 17.35—20.93 
24.66 一 28.90 Wm 20.93 一 25.77 


N 
入 


40 一 60cm ~ 集 水 线 
土壤 水 分 含量 /% 坡度 
mm 5.661—7.272 <15° 
BN 7.272—8.645 15° 25° 
Ba 8.645 一 9.815 口 >25° 
mn 9.815—10.81 坡 向 
ml 10.81—11.98 国 阴 坡 
Ba 11.98—13.35 襄阳 坡 
a 13.35 一 14.97 海拔 /km 
a 14.97—16.86 1.9—2.1 
ma 16.36—19.07 2.1 一 2 2 
ma 19.07 一 21.68 国 2.2 一 2 3 
国 2.3 一 2 .5 


图 3 不 同 土壤 深度 土壤 水 分 的 空间 分 布 特征 及 地 形 因子 又 加 
Fig.3 Spatial distribution characteristics of soil moisture content across depths and overlay of topographic factors 


SMC (80il] moisture content) 表示 土壤 水 分 含量 


从 表 5 中 可 以 看 出 ,不 同 次 度 的 土壤 水 分 含量 与 地 形 因子 (海拔 .坡度 . 坡 向 ) 表 现 出 显著 的 相关 关系 (P< 
0.01 ) ,而 与 植被 指数 NDVI 的 负 相 关 关 系 , 却 不 显著 ,说 明山 区 森林 土壤 水 分 的 空间 异 质 性 主要 是 由 地 形 因 
素 引 起 的 。 坡 向 指数 与 去 层 土壤 水 分 随 土 层 表现 出 逐渐 增强 的 负 相 关 关 系 , 说 明生 境 为 阴 坡 、 湿 冷 的 环境 , 蒜 
发 量 较 小 ,有 利于 土壤 的 蕃 水 。 海 拔 对 0 一 60cm 土壤 水 分 的 影响 为 强烈 的 正 相 关 性 , 即 海拔 越 高 ,土壤 水 分 
越 易 保持 。 坡 度 越 大 ,土壤 含 水 量 越 低 ,反之 , 越 高 。 海 拔 和 坡度 对 土壤 水 分 的 影响 随 着 土 层 的 加 深 逐 渐 减 弱 
(海拔 :0.73 亿 0.69) 坡 度 :-0.58 一 -0.50,P<0.01) 。 另 外 ,土壤 水 分 与 土壤 化 学 性 质 ( 有 机 碳 、 全 所 含量 ) 具 有 
显著 的 正 相 关 关 系 (P<0.01) ,说 明 土壤 有 机 质 和 全 所 含量 对 土壤 水 分 具有 显著 影响 。 结 合 植被 分 布 情况 , 土 
壤 水 分 的 高 值 斑 块 区 域 主要 分 布 在 华北 落叶 松 近 熟 林 下 ,而 落叶 松 成 熟 林 、 混 交 林 和 灌木 林 下 土壤 水 分 值 相 
对 较 低 。 


表 5 不 同 深度 土壤 水 分 与 地 形 、 植 被 因子 和 土壤 化 学 性 质 之 间 的 相关 关系 


Table S Correlations between soil moisture content, topographic, vegetation factors and soil chemical properties 


土 层 坡 向 指数 海拔 坡度 植被 指数 。 土壤 有 机 碳 仿 量 土壤 全 握 含 量 

Depth /em TRASP Altitude/m Slope/ (°) NDVI Soil organic carbon/ % Soil totalnitrogen/ % 
0 一 20 一 0.406 0.733 一 0985 -0.156 0.601 一 0.685 

20 一 40 一 0.413 0.711 一 0.380 一 -0.210 0.571 一 0.649 一 

40 一 60 一 0.5305 0.695 一 一 0.504 一 一 0.227 0.671 一 0.762 一 


x** P<0.01 水 平 上 显著 相关 
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3.1 小 流域 土壤 水 分 的 时 间 稳 定性 与 变异 系数 

有 学 者 提出 在 研究 土壤 水 分 的 空间 分 布 时 , 因 土 壤 含 水 量 在 很 大 程度 上 取决 于 土壤 结构 ,假设 其 在 时 间 
上 具有 一 定 的 稳定 性 , 且 用 非 参数 法 和 时 间 标 准 均 方差 比较 进行 了 验证 ,证 明 土 壤 水 分 的 空间 分 布 具有 较 强 
的 时 间 稳 定性 5 。 本 研究 取样 的 间隔 时 间 较 短 , 天 气 状况 稳定 ,因此 土壤 水 分 具有 相对 的 时 间 稳 定性 。 暖 温 
带 亚 高 山 森 林 小 流域 土壤 水 分 各 层 变异 系数 介 于 32.62% 一 38.98% 之 间 , 此 值 范围 稍 高 于 阔 叶 红 松林 《各 % 碟 
右 ) 所 及 人 工 油 松林 (10% 一 25% ) 5 ,可 能 是 由 于 已 有 的 山地 森林 土壤 水 分 空间 异 质 性 研究 大 多 基于 样 方 和 
坡 面 尺度 ,森林 小 流域 尺度 的 范围 较 大 ,导致 植被 .土壤 及 地 形 类 型 有 所 差异 ;另外 砾石 含量 相对 较 多 且 采 样 
节点 土壤 较 干 (12.54% 一 15.75% ) 也 是 引起 变异 系数 稍 大 的 原因 之 一 。 土 壤 水 分 变异 系数 的 层次 递减 规律 与 
莫 菲 等 2 沁 对 相近 纬度 的 六 盘山 洪 沟 小 流域 落叶 松 植被 下 土壤 水 分 的 研究 结果 相 侧 。 
3.2 ”小 流域 土壤 水 分 的 空间 自 相 关 性 

研究 区 土壤 水 分 表现 为 强烈 的 空间 自 相 关 性 ,其 结构 性 方差 在 81.40% 一 91.30%, 变 程 在 1.1 一 1.4 km 之 
间 。 和 森林 土壤 水 分 在 时 间 稳 定 下 较 强 的 空间 自 相关 性 与 已 有 的 多 数 研究 结果 相 只 合 "“”]。 但 Western 等 提 
出 当 取 样 间 隔 大 于 变 程 .取样 点 数量 太 少 以 及 土壤 水 分 的 测定 误差 居于 其 变异 等 因素 时 会 引起 土壤 水 分 弱 或 
无 空间 相关 性 :2 。 与 已 有 研究 相 比 ,本 文 的 最 大 变 程 相对 小 流域 尺度 的 其 他 研究 要 大 '2 ,分 析 其 原因 一 方 
面 可 能 是 接 石 沟 小 流域 分 布 着 大 面积 的 华北 落叶 松 ( 覆 盖 面 积 约 占 80%) ;树木 的 年 龄 .空间 分 布 和 垂直 结构 
都 相对 均一 ,植被 斑 块 均 质 化 ,使 得 土壤 水 分 的 自 相 关 范 围 延 恬 。 通 可 能 也 与 地 形 特征 有 关 , 研 究 区 范围 的 平 
均 坡 长 在 900 m 左右 ,在 一 定 程度 上 可 以 解释 空间 自 相 关 变 化 的 范围 。 值 得 注意 的 是 , 现 有 的 采样 尺度 可 能 
在 更 大 范围 内 仍然 存在 空间 异 质 性 。 本 研究 的 平均 采样 距离 为 263 m, 低 于 变 程 范围 的 一 半 , 说 明 在 此 小 流 
域 尺 度 下 ,采样 密度 和 方法 足以 说 明 该 区 域 土壤 水 分 的 空间 异 质 性 ,该 结果 对 今后 典型 森林 小 流域 的 取样 间 
隔 设 置 具有 一 定 的 借鉴 作用 。 
3.3 不 同 层次 土壤 水 分 空间 异 质 性 的 影响 因素 

土壤 水 分 空间 异 质 性 与 植被 .地形 因子 和 土壤 性 质 有 着 密切 的 联系 。 研 究 结 果 得 出 ,土壤 含水 率 随 着 海 
拨 的 升 高 显著 上 升 。 结 合 植被 分 布 诲 拔高 的 地 区 植被 多 为 华北 落叶 松 人 工 纯 林 ,而 海拔 相对 低 的 地 区 还 有 
一 些 阔 叶 林 分 布 ,如 白桦 ( Betwla Platyphylla) ,黄花 柳 ( Salix caprea) 等 ,阔叶树 种 的 叶片 表面 积 大 ,蒸发 量 较 
大 ,导致 土壤 水 分 偏 低 和 而 海拔 高 的 针 叶 林 可 以 减少 水 分 和 热量 的 蒸发 ,因而 土壤 水 分 偏 高 。 男 外 ,坡度 越 陡 ， 
人 渗 率 低 ,下 渗 和 径流 越 强 ,导致 土 壤 水 分 含量 较 少 ,反之 越 易 于 保持 水 分 。 坡 向 偏向 于 阴 坡 , 落 发 量 小 ,有 利 
于 水 分 聚集 ,反之 藻 发 量 大 ,土壤 水 分 损耗 大 ,与 一 般 研究 的 结论 相似 ”1 。 在 与 植被 的 关系 方面 , Schume 
等 ' 引 认为 森林 强大 深厚 的 根系 对 于 土壤 水 分 的 吸收 作用 比 农田 或 草地 等 生态 系统 更 明显 ,是 影响 土壤 水 分 
的 主 控 因 了 予 % 而 本 人 研究 得 出 植被 NDVI 指数 对 于 土壤 水 分 有 弱 的 负 影 响 。NDVI 意味 着 高 的 叶 面 积 指数 和 植 
被 覆盖 谋 , 相 应 的 获 发 量 也 越 大 ,导致 与 土壤 水 分 负 相 关 , 且 在 采样 时 间 节 点 ,土壤 较 干 ,地 形 和 土壤 性 质 的 影 
响 要 为 于 植被 ,这 与 Grayson 等 小 的 研究 结论 一 致 。 从 具体 林 分 结构 看 , 近 熟 林 的 冠 层 不 如 成 熟 林 的 郁 闭 , 冠 
幅 、 冠 高 和 根系 深度 均 小 于 成 熟 林 , 叶 片面 积 又 较 混 交 林 的 要 小 ,导致 近 熟 林地 的 蒸发 蒸腾 小 于 成 熟 林 与 混交 
林 而 灌木 林 分 布 分 散 ,其 所 在 区 域 土 层 较 浅 , 不 易 储存 水 分 ,说 明 林 分 的 结构 特征 也 是 导致 土壤 水 分 空间 异 
质 性 的 重要 原因 。 另 外 ,有 研究 表示 土壤 的 持 水 性 能 还 与 有 机 质 含量 及 大 于 0.25 mm 的 水 稳 性 团聚 体 有 
关 !21 ,这 一 观点 与 本 文 森林 土壤 水 分 含量 和 有 机 碳 、 全 所 显著 正 相 关 的 结果 相 吻 合 。 从 土壤 不 同 层 次 看 ， 
Rosenbaum 等 中 研究 表明 ,土壤 上 层 (<20cm) 水 分 主要 受到 坡 位 坡度 等 地 形 条 件 和 小 气候 条 件 下 的 蒸发 力 
影响 , 随 着 士 层 加 深 , 土壤 水 分 受 莹 发 作用 减 小 ,而 受 植被 蒸腾 作用 加 大 。 本 文 统 计 结 果 显 示 ,海拔 和 坡度 等 
地 形 因 子 在 随 着 土 层 深度 加 深 , 与 土壤 水 分 的 相关 系数 减弱 ,而 植被 指数 (NDVI) 则 与 土壤 水 分 的 相关 性 增 
强 , 由 此 可 看 出 , 随 着 土 层 的 加 深 , 地形 的 作用 逐渐 减弱 ,植被 蒸腾 作用 加 大 。 
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在 亚 高 山 山地 华北 落叶 松林 分 布 为 主 的 小 流域 ,地 形 较 复杂 ,土壤 水 分 含量 及 变异 系数 在 表层 (0 一 
20cm) 最 大 , 随 着 土屋 加 深 逐 渐 递 减 。 三 层 深度 土壤 水 分 均 呈 强烈 的 空间 自 相 关 , 其 中 以 中 间 层 的 空间 自 相 
关 性 最 强 , 由 结构 因素 部 分 引起 的 空间 自 相 关 性 占 主导 地 位 ,随机 部 分 引起 的 误差 较 小 。 采 样 密度 能 满足 土 
坏 水 分 空间 变异 分 析 的 要 求 。 在 植被 覆盖 异 质 性 小 的 山地 ,地 形 起 伏 对 土壤 水 分 的 影响 起 着 主导 作用 ,但 随 
着 土 层 的 加 深 , 其 作用 减弱 ,植被 特征 ( 林 分 年 龄 与 结构 ) 的 影响 逐渐 增强 。 土 壤 水 分 含量 随 着 海拔 的 呈 升 
坡度 的 减 小 和 坡 向 的 偏 北 而 增高 ,反之 亦 然 。 地 形 的 影响 与 植被 因素 密 不 可 分 ,土壤 水 分 从 山顶 华北 落叶 松 
近 熟 林 到 沟 口 混交 林 、 灌 木林 呈 梯 度 递减 趋势 。 

研究 土壤 水 分 空间 异 质 性 的 最 大 困难 在 于 土壤 水 分 随时 间 变 化 较 大 , 且 在 不 同 尺 度 上 变化 规 带 和 不同” 。 
本 文 只 对 特定 时 间 内 土壤 水 分 空间 变异 性 的 一 些 规律 进行 探讨 ,要 全 面 了 解 华北 地 区 典型 山地 和 木林 土壤 水 分 
的 空间 变异 规律 ,还 需要 对 土壤 水 分 随时 间 变 化 进行 观测 ,并 在 不 同 尺度 上 做 进一步 研究 。 
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